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We present a new method to determine the momentum dependence of the N to A transition form factors and 
demonstrate its effectiveness in the quenched theory at (3 = 6.0 on a 32 3 x 64 lattice. We address a number of 
technical issues such as the optimal combination of matrix elements and the simultaneous overconstrained analysis 
of all lattice vector momenta contributing to a given momentum transfer squared, Q 2 . 



1. Introduction 

The N to A transition form factors encode 
important information on hadron deformation 
and have been studied carefully in recent experi- 
ments 1 . In this work we present the first lattice 
evaluation of the momentum dependence of the 
magnetic dipole, Ml, the electric quadrupole, E2, 
and the Coulomb quadrupole, C2, transition am- 
plitudes. They are calculated in the quenched ap- 
proximation on a lattice of size 32 3 x 64 at (3 = 6.0 
with Wilson fermions with sufficient accuracy to 
exclude a zero value of E2 and C2 at low Q 2 . 
This accuracy is achieved by applying two novel 
methods: 1) We use an interpolating field for the 
A that allows a maximum number of statistically 
distinct lattice measurements contributing to a 
given Q 2 . 2) We extract the transition form fac- 
tors by performing an overconstrained analysis of 
the lattice measurements using all lattice momen- 
tum vectors contributing to a given Q 2 value 0- 

2. Evaluation of the three-point function 

The evaluation of the three-point function 
QAj N (t 2 ^t 1 ]p\p]T) can be done either using 
the fixed current approach as in previous lat- 
tice calculations |3I4) or the fixed sink approach, 
where the current can couple to the backward 
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sequential propagator at any time slice t\ car- 
rying any lattice momentum allowing the 
evaluation of the form factors at all possible 
momentum transfers. As shown in Fig. ^ for 
the same statistics, the errors in Ml in the 
fixed sink method are almost three times as 
large as those obtained in the fixed current ap- 
proach. Therefore in order to make use of the 
advantages of the fixed sink approach, we must 
first reduce the errors. We start by modify- 
ing the ratio used in ref. 0] so that we do not 
need to evaluate both {G^ N (t2, fa; p ', p; F)} 
and (G^i A (t2,fa;p ', p; T)}, since this would re- 
quire two inversions. Instead we use the ratio 
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in the notation of ref. 0]. We use kinematics 
where the A is produced at rest and so q = 
p' p = —p. We &xt2/a= 12 and search for a 
plateau of R a (t2, fa; p ', p ; T; fx) as function of fa. 

In the fixed sink approach, the index a of the 
A and projection matrix T are fixed and therefore 
we need to determine the most suitable choice of 
three-point functions from which to extract the 
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Figure 1. Gmi in Bohr magnetons at n = 0.1554. 

Sachs form factors Gm\,Ge2 and Gc2 0]- For 
example Gmi can be extracted from 

n ff ( q ; r 4 ;/i) = l v4 £ ff4 wyeAfi(Q 2 ) (2) 



where A is a kinematical coefficient and T = 
r 4 = \ i q M in Eq. JTJ. This leaves 3 

choices for er i.e there are three statistically inde- 
pendent matrix elements yielding Gmi, each re- 
quiring a sequential inversion. However, due to 
the e factor, fixing a means that only momentum 
transfers in the other two directions contribute. 
Instead, if we take the symmetric combination, 
Si(q;^) = X)tr=i n CT (q;r 4 ;/i), momentum vec- 
tors in all directions contribute. This combina- 
tion, which we refer to as sink S\ , is built into the 
A interpolating field and requires only one inver- 
sion. To take full advantage of the number of lat- 
tice vectors contributing to a given Q 2 we perform 
an overconstrained fit by solving the overcomplete 
set of equations P(q; /i) = -D(q; /tx) • F(Q 2 ) where 
P(q; /i) are the lattice measurements of the ra- 
/ Gmi \ 

tio of Eq. P, F = Ge2 and, with N be- 

ing the number of current directions and momen- 
tum vectors contributing to a given Q 2 , D is an 
N x 3 matrix which depends on kinematical fac- 
tors. We extract the form factors by minimizing 

X 2 = Ef=i ^ (E]=i DkjFj - PfeJ , where w k 
are the errors in the lattice measurements, us- 
ing singular value decomposition of D. In Fig.^ 
we compare the results for Gmi using an over- 
constrained analysis with sink type Si to our old 
analysis with fixed a = 2 and /i = 3. The errors 
with our new analysis are reduced at all values of 
Q 2 and are now equal to the error obtained using 
the fixed current approach. 
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Figure 2. Top: Gmi, middle: Ge2, bottom: 
Gc2 in Bohr magnetons at K = 0.1554 using 
n3(q;ri;j) (crosses), n 3 (q;ri;j) & ni(q;r 3 ;j) 
(filled triangles), Si or 5*2 (asterisks) and S3 
(open circles) and an overconstrained analysis. 

Using = i { a ^ q ^ instead of T 4 gives 

another six statistically independent 3-point func- 
tions from which Gmi can be extracted: 

U a (q;T k ;j) = b\ i {p a 5 kj ■. - p k S aj ) Gmi(Q 2 ) 
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where B and C are kinematical coefficients. To 
isolate the benefits of using S\ we compare in 
Fig. El Gmi obtained using Si to the the ones ob- 
tained by fixing a = 3 and Y\ in Eq. © 0]. All 
results are now obtained with the overconstrained 
analysis using 50 configurations. As can be seen, 
Si produces results with the smallest errors and 
it is therefore the optimal sink for Gmi- 

For the extraction of the quadrupole mo- 
ments, we consider the symmetric combination 
S 2 (q; n) = J2ajik=i ^(q; r fe ; fi) from which both 
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Qe2 and Qc2 can be extracted when the current 
is in the spatial direction. When the current is 
in the time direction, S2 provides a statistically 
independent way for evaluating Gc2, at no ex- 
tra cost. Another combination to extract the 
quadrupole form factors is S3 = n 3 (q ;T 3 ;/j,) — 
|(lli(q +II 2 (q ; T 2 ; n)) , which, unlike 52, 

contributes at the lowest value of Q 2 . As can be 
seen in Fig. [21 £2 produces results with smaller 
errors as compared to those using H^(Ti;j) and 
TJi(Ts;j) for both E2 and C2. Increasing the 
the statistics from 50 to 200 configurations brings 
agreement between S2 and S3 for Gc2 as well. 

3. Results and Conclusions 

We analyse 200 configurations at re = 0.1554, 
0.1558 and 0.1562 corresponding to m^/nip = 
0.64, 0.59 and 0.50 respectively. We use the nu- 
cleon mass at the chiral limit to set the lattice 
spacing a, obtaining a" 1 = 2.04(2) GeV. Using 
the optimal sink Si our results for G%u are shown 
in Fig. as a function of Q 2 , where 
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Results in the chiral limit are obtained by per- 
forming a linear extrapolation in m 2 .. On 
the same figure, we also show the exper- 
imental values as extracted from the mea- 
sured cross sections using the phenomenolog- 
ical model MAID 0. Although the lattice 
data in the chiral limit lie higher than the 
MAID data both data sets are well described 
by the the phenomenological parametrization 

G a {Q 2 ) = Ga{0)(l + aQ 2 )eM-lQ 2 )G p E {Q 2 ) 
for a = Ml, £2 and C2 and G P E (Q 2 ) = 1/(1 + 
Q 2 /0.71) 2 is the proton electric form factor. 
These fits are shown in Fig. [31 by the solid lines. 
In Fig. [3 we show the ratios EMR = Rem = 



GE2(q 2 ) 
<?Afl(? 2 



and CMR 



d _ q| Qcijq ) 
n SM - 2mA g M1 ( q 2)- 



Results in the chiral limit are obtained by per- 
forming a linear extrapolation in m 2 . As ex- 
pected, both EMR and CMR become more neg- 
ative as we approach the chiral limit. The re- 
sults for EMR are in agreement with experimental 
measurements whereas CMR is not as negative as 
experiment at low Q 2 . We believe that CMR is 



particularly sensitive to the absence of sea quarks 
and thus a good probe of unquenching effects. 
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Figure 3. Top 
CMR as function of Q 2 at n 
k = 0.1558 (open triangles), 



0.6 o.f 
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Gmi, middle: EMR, bottom: 
= 0.1554 (crosses), 
k = 0.1562 (aster- 



isks) and in the chiral limit (filled circles). Filled 
triangles show G\ji extracted from measurements 
using MAID jfjj. The dashed line is a fit to the 
lattice data using aexp(— bQ 2 ). 
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